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ABSTRACT 

A  fixed  base  helicopter  simulation  was  implemented 
which  may  be  easily  altered  to  represent  a  variety  of  con- 
ventional single  rotor  aircraft.   Test  subjects  performed 
90-second  tracking  tasks  utilizing  two  cockpit  displays  for 
longitudinal  control  during  landing  approach.   The  performance 
data  generated  were  compared  with  predicted  values  from  an 
optimal  pilot  modeling  technique.   Tentative  conclusions 
were  drawn  as  to  the  effectiveness  of  each  display  as  well 
as  the  predictive  power  of  the  pilot  model. 
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g    Acceleration  of  gravity,  ft/sec  . 

I    Aircraft  mass  moment  of  inertia  about  y  stability  axis, 

2 

slug-ft  , 

M  Moment  about  y  stability  axis,  ft-lbs. 

m  Aircraft  mass,  slugs. 

q  Aircraft  angular  rate  about  y  stability  axis,  radians/sec. 

U  Aircraft  reference  flight  velocity,  ft/sec, 

u    Perturbation  flight  velocity  along  x  stability  axis,  ft/sec 
u    Horizontal  wind  gust  velocity,  ft/sec, 

w    Perturbation  flight  velocity  along  z  stability  axis, 
ft/sec. 

w    Vertical  wind  gust  velocity,  ft/sec. 

X    Force  component  along  x  stability  axis,  lbs. 

Z    Force  component  along  z  stability  axis,  lbs. 

6,    Cyclic  pitch  control  input,  ft.  displacement  at  the 
pilot's  hand. 

6    Collective  pitch  control  input,  ft.  displacement  at 
the  pilot's  hand. 

9     Perturbation  pitch  angle  about  y  axis,  radians. 

The  following  stability  derivatives  are  defined  for  a 
reference  flight  condition  in  stability  axes. 
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I.   INTRODUCTION 

The  use  of  vertical  take-off  and  landing  aircraft  has 
increased  significantly  in  both  military  and  civil  applica- 
tion.  Helicopter  flight  under  instrument  flight  rules  has 
been  accomplished  on  a  routine  basis  for  quite  a  while  but 
low  speed  and  terminal  operations  still  require  visual 
references  outside  the  aircraft.   The  military  services  have 
had  some  success  in  low  speed  instrument  flight  but  it  is 
generally  conceded  that  much  work  needs  to  be  done  in  the 
area  of  cockpit  information  systems.   One  popular  approach 


has  been  the  use  of  electronic  displays  utilizing  cathode- 
ray-tubes  which  have  the  advantage  of  almost  unlimited 
flexibility  in  format  and  symbolism  (References  1  and  2). 
Additionally,  extensive  work  has  been  done  in  the  field  of 
pilot  modeling  and  pilot-in-the-loop  system  performance  (Ref.3); 
but  in  order  to  optimize  pilot  performance  it  is  necessary 
to  consider  the  nature  and  extent  of  information  presented 
to  the  pilot.  \  This  project  was  aimed  at  the  construction 
of  a  fixed  base  helicopter  flight  simulator  and  the  comparison 
of  real  pilot  performance  with  pilot  model  predictions  under 
a  limited  range  of  conditions  and  using  two  visual  displays. 
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II,   METHOD  OF  INVESTIGATION 

A  helicopter  simulation  was  constructed  to  represent 
the  longitudinal  flight  dynamics  of  an  unaugmented  (i.e., 
without  the  usual  mechanical  stabilization  bar)  UH-1H 
helicopter.   A  fore  and  aft  moving  control  stick  provided 
cyclic  pitch  control  inputs  and  a  collective  pitch  level 
produced  power  commands.   A  hybrid  computer  gave  real  time 
solutions  to  the  helicopter  equations  of  motion,  generated 
each  of  two  cockpit  displays,  and  computed  performance 
data.   Performance  values  for  each  display  mode  were 
analytically  predicted  using  optimal  pilot  modeling  tech- 
niques and  compared  with  simulation  data  for  two  test 
subjects. 

III.   APPARATUS  AND  PROCEDURES 

The  longitudinal  motion  of  a  helicopter  can  be  described 
by  the  following  equations  of  motion: 


u  =   Xu   +   Xw  +   Xq-gO+X~      6.     +  X.         6      -Xu      -  X  w 
u  w  q^  6,       b  6  c  u  g  w  g 


w=Zu   +    Zw+(U      +Z)q   +    Z.      6.     +   Zr      6      -Zu      -Zw 
u  w  o  q   ^  o        b  6        c  ug  wg 


q   =    (M   +M.Z    )u   +    (M   +M.Z    )w   +    (M   +M.(U    +Z    ))q    +    (M.    H-M.Z^     )6,     + 
vuwu  v    w     w  w'  vqwvo      q  5b  b 

+    (Mx    +M.Z-    )6      -    (M   +M.Z    )u      -    (M   +M.Z    )w 
o.      w   6    '    c  u      wuyg        vw     ww'g 
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0  =  q 

h  =  -w  +  U  0 
o 

6,  represents  cyclic  control  input,  6   collective  input, 

and  u   and  w  horizontal  and  vertical  wind  gust  velocities, 
g      g 

Table  I  lists  the  stability  derivative  values  used  in  this 
simulation.  The  equations  above  incorporate  the  following 
simplifying  assumptions: 

1.  The  vehicle  has  been  idealized  as  a  rigid  airframe  to 
which  is  attached  a  rotor. 

2.  The  rotor  is  described  by  its  tip  path  plane  whose 
orientation  determines  propulsive  and  aerodynamic  forces 
and  moments. 

3.  No  rotor  degree  of  freedom  is  considered  other  than 
control  inputs  which  serve  to  describe  instantaneous 
tip  path  plane  orientation. 

4.  Coupling  between  longitudinal  and  lateral  motion  has 
been  ignored. 

5.  Motion  has  been  linearized  as  small  perturbations  about 
a  predetermined  reference  flight  path, 

The  hybrid  computer  used  consisted  of  a  Scientific  Data 
Systems  SDS  9300  digital  computer,  A  Comcor  CI-5000  analog 
computer,  and  an  Adage  AGT-10  stroke  written  graphics 
processor,  as  illustrated  in  Figure  1.   The  digital  computer 
was  used  to  control  the  analog  and  graphics  processor  and 
to  generate  performance  data,   A  program  listing  is  found  in 
Appendix  A.   In  addition  to  the  two  display  modes  to  be 


12 


described,  provision  has  been  made  for  the  inclusion  of 
flight  director  symbolism  in  a  third  display  mode  for  future 
study.   The  analog  computer  gave  solutions  to  the  equations 
of  motion  in  real  time  (Appendix  B) ,  and  the  graphics  processor 
displayed  textual  and  graphic  information  on  a  cathode  ray 
tube  in  front  of  the  test  subject.   Pilot  control  inputs 
were  made  directly  to  the  analog  using  electrical  voltage 
signals.   The  spring  restrained  cyclic  stick  operated  a  gear 
driven  potentiometer  and  was  linear  in  both  force  and  dis- 
placement.  The  collective  lever  was  unrestrained  except  for 
a  friction  brake  and  also  used  a  gear  driven  potentiometer 
to  generate  a  signal  linear  in  displacement.   The  friction 
brake  required  a  slight  breakout  force.   Figure  2  shows  the 
physical  arrangement  of  the  apparatus. 

Atmospheric  turbulence  was  modeled  by  a  sum  of  sine 
waves : 

5 

u   =   E  A.sinoo.t,  where     A.,  =  4.472         oj-  =  0.070 
g    i=1  i     i  1  1 

A2  =  3.536  w2  =  0.210 

A3  =  2,236  0)3  =  0.350 

k  =  2,738  co   =  0.768 

Ac  =  2.236  ojc  =  1.765 

and  similarly, 

5 

w     =      E   A.sinoo.t,    where  A-    =   4,472  oj1    =   0.420 

g        i=1    i  i  l  l 

A2   =    3.536  u>2   =    1.260 

A3   =   2.236  oj3   =    2.100 

A.    =    2.738  oj„    =   4.608 

4  <\ 

A5   =    2.236  oo5   =    10.050 
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The  amplitudes  and  frequencies  above  were  chosen  so  that 
the  integration  over  oo  of  the  power  spectral  densities 
closely  approximated  that  for  the  turbulence  spectra  shown 
in  Table  II  (Ref,  4).   The  advantage  of  this  turbulence 
representation  is  that  it  is  repeatable,  has  a  mean  value 
of  zero,  and  appears  completely  random  to  the  test  subject. 

Two  different  displays  were  presented  to  the  test 
subjects  for  evaluation.   The  first  was  basically  a  real 
world  display  similar  to  current  instrument  installations  and 
consisted  of  a  horizon  bar  and  aircraft  reference  for  pitch 
control  along  with  a  glideslope  deviation  indicator.   Air- 
speed deviation  was  provided  by  angular  displacement  of  the 
aircraft  reference  "wings."   The  second  display  had  symbolism 
identical  to  the  baseline  display  but  was  quickened  with  rate 
information  according  to  the  equations: 

0   .  .  =  0  +  k  »0 
quick        8 

h   .  ,  =  h  +  k«h   . 
quick        h 

The  quickening  gains,  k»  and  k • ,  were  simple  pure  gain 

approximations  to  the  pilot  transfer  functions  obtained  in 
the  modeling  procedure.   The  display,  shown  in  Figure  3,  was 
6  1/2  inches  wide  and  7  inches  high.   The  nominal  eye-to- 
display  distance  was  2  1/2  feet. 

After  considerable  informal  training,  each  subject 
received  a  formal  training  session  in  which  he  used  the  base- 
line and  quickened  display  for  90  second  runs  each.   In  the 
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formal  data  sessions,  each  subject  used  each  display  for  ten 
90  second  runs.   The  data  sessions  went  as  follows: 

The  baseline  was  flown  for  seven  runs,  followed  by  the 
quickened  for  seven.   Then  the  quickened  was  flown  for 
three  runs  followed  by  the  baseline  for  three.   Five  of 
the  "best"  runs  for  each  display  were  then  selected.   For 
ten  data,  runs  with  each  display,  mean  scores  were  computed 
for  each  of  the  root  mean  square  performance  measures 
(u,  9,  9,  h,  <5,  ,  6  ).   Then  each  RMS  score  for  each  display 
and  run  was  normalized  by  dividing  it  with  the  respective 
mean  RMS  score.   The  normalized  RMS  scores  were  then  summed 
to  provide  a  single  scalar  index  of  performance  for  each 
display  and  run. 
The  data  in  Table  IV  represent  the  means  and  standard 
deviations  of  the  RMS  scores  for  the  runs  with  the  five 
lowest  indices  of  performance  for  each  display. 

Figures  5  and  6  graphically  portray  the  more  pertinent 
performance  data  of  Table  IV  for  subject  A  (non-pilot)  and 
subject  B  (pilot).   Also  indicated  are  the  model  predictions 
which  best  correlate  with  the  data,  i.e.  F  =  0,25. 

IV.   OPTIMAL  PILOT  MODEL  AND  PREDICTIONS 

The  analytic  pilot  model  is  shown  schematically  in 
Figure  4,  and  is  basically  the  same  as  that  offered  by 
Kleinman  and  Baron  (Ref.  5).   This  optimal  pilot  model 
includes  the  following  assumptions: 
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1.  Subject  to  inherent  limitations,  the  well  trained  and 
motivated  pilot  behaves  in  an  optimal  manner. 

2.  With  certain  modifications,  the  solution  to  an  optimal 
linear  control  and  estimation  problem  will  predict  the 
pilot's  control  characteristics.   These  modifications 
are : 

a.  A  pure  time  delay  is  included  in  each  of  the  pilot's 
control  outputs. 

b.  Each  neuromuscular  output  system  may  be  modeled  as 
a  first  order  lag. 

c.  Each  variable  that  the  pilot  observes  is  assumed  to 
contain  pilot  induced  additive  white  noise  which 
scales  with  the  variance  of  the  observed  variable. 

d.  Each  control  output  is  assumed  to  contain  pilot 
induced  additive  white  noise  which  scales  with  the 
variance  of  the  control. 

e.  The  pilot  is  able  to  perceive  the  first  derivative 
of  a  variable,  also  noise  contaminated,  but  no 
higher  derivative. 

The  index  of  performance  for  the  optimization  procedure  was 

chosen  subjectively  by  the  analyst  to  represent  the  task 

objectives  as  perceived  by  the  pilot.   Model  performance 

predictions  are  given  in  Table  III  for  each  of  the  displays. 

The  parameter  F  represents  the  fraction  of  the  pilot's 

attention  devoted  to  the  control  task  as  a  whole.   For  example, 

F  =  0.25  suggests  that  25%  of  the  pilot's  attention  is  devoted 

to  control  and  75%  is  dedicated  to  non-control  activity  such 

as  communication,  etc, 
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V.   RESULTS 

Table  IV  lists  the  root  mean  square  performance  data. 
Figures  5  and  6  give  a  graphical  comparison  of  the  most 
pertinent  performance  measures  and  the  model  predictions  that 
best  correlate  with  the  data  (F  =  0.25).   For  each  subject, 
the  mean  value  is  indicated,  as  well  as  plus  and  minus  one 
standard  deviation.   In  each  case  the  pilot  model  predicted 
performance  with  a  fair  degree  of  accuracy.   Additionally, 
the  model  predicted  improvement  for  each  performance  measure 
using  the  quickened  display,  and  the  test  subjects  verified 
this  prediction.   However,  the  magnitude  of  performance 
improvement  was  greater  in  each  case  than  was  anticipated. 
It  is  possible  that  the  task  as  perceived  by  the  pilot  may 
not  be  the  same  for  each  display  mode.   The  task  as  understood 
may  be  to  keep  the  deviation  symbols  within  a  finite  distance, 
say  1/8",  of  their  nominal  position,  and  indeed  there  may 
exist  a  threshold  of  error  below  which  the  pilot  makes  no 
control  input  at  all.   Since  the  gains  driving  the  deviation 
symbolism  vary  considerably  between  modes,  each  display  then 
defines  a  new  task.   One  possible  arrangement  would  be  to 
provide  additional  symbolism  for  the  task  variables,  i.e. 
airspeed  and  glideslope  deviations,  which  would  be  common 
between  modes. 

The  dramatic  reduction  in  cyclic  motion,  o,  ,  was  not 
predicted  by  the  model.   Typical  time  histories  of  cyclic 
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motion,  Figures  7  and  8,  show  the  nature  of  this  improvement. 
Again  the  task  defined  by  the  display  mode,  in  this  case 
pitch  attitude  control,  may  be  different.   For  both  subjects, 
glideslope  deviation  control  was  poorest  in  comparison  with 
model  predictions,  which  suggests  that  more  pilot  attention 
was  given  to  attitude/airspeed  control  than  to  the  glideslope. 
Furthermore,  collective  motion  was  less  for  both  subjects 
than  was  predicted.   While  the  cyclic  input  was  of  a  continuous 
nature,  collective  inputs  were  a  series  of  steps  (Figures  7 
and  8)  in  each  case.   Such  control  action  probably  represents 
a  wait  and  see  technique  used  by  pilots  as  a  result  of  the 
lag  in  collective  response.   This  phenomenon  may  in  part 
explain  the  discrepancy  between  predicted  and  measured 
glideslope  deviations  since  the  optimal  control  model  did  not 
address  this  type  of  control  motion. 

Pilot  comments  indicated  that  the  quickened  display 
was  more  fatiguing  than  the  baseline  and  represented  an 
increase  in  pilot  workload.   Unfortunately,  no  measure  of 
pilot  workload,  such  as  number  of  control  inputs  per  minute, 
was  recorded.   Furthermore,  no  objective  measure  of  the  "well 
trained  and  motivated"  assumption  was  available  for  each 
subject . 

VI,   CONCLUSIONS 

The  use  of  quickening  in  the  cockpit  display  significantly 
reduced  control  motion  and  improved  pilot  performance.   The 
predictive  power  of  the  optimal  pilot  model  seemed  adequate 
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even  though  this  report  represents  only  a  first  iteration. 
The  quantitative  quality  of  the  model  should  improve  as 
adjustments  are  made.   The  model  index  of  performance  can 
be  adjusted  to  account  for  the  pilot's  apparent  greater 
concern  with  pitch  attitude/rate  deviations  as  compared  with 
altitude/rate  deviations.   In  addition,  the  displays  could 
be  revised  to  eliminate  any  variation  in  the  task  presented 
to  the  pilot.   As  the  control  theoretic  pilot  model  is 
refined  and  validated  through  testing  it  should  become  a 
valuable  tool  in  the  development  of  cockpit  information 
systems  and  the  expanded  utility  of  the  helicopter. 
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TABLE  I 


UH-1H  LONGITUDINAL  DERIVATIVES  AS  USED  IN  SIMULATION 
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Z 

w 

Z 

q 

M 
u 

M 
w 

M 

q 

M. 
w 

X6b 

6c 

Z6b 

Z6c 

M6b 

M6c 


-0.11985 

1/sec 

0.35800 

1/sec 

1.9359 

ft/ sec 

0.28007 

1/sec 

-1.0775 

1/sec 

-1.4597 

ft/sec 

0.0058048 

1/sec-f t 

-0.0030342 

1/sec-ft 

-0.45465 

1/sec 

0.0 

1/ft 

-3.4680 

2 
1/sec 

59,688 

2 
1/sec 

48.486 

1  .    2 
1/sec 

-176.13 

2 
1/sec 

-1,5757 

2 

1/ft-sec 

0.0 

2 
1/ft-sec 
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TABLE  II 
TURBULENCE  MODELS 


Vertical  Turbulence 

$       (CO) 

w  w  v  ' 

g  g 

2 

2awLw     1       ^2    4#    2 
=  — 7j —  = ft   rad/sec 

1 +(tr^) 

o 

a 
w 

=  5  ft/sec 

L 
w 

=  100  ft 

U 
o 

=  60  knots  (101.3  ft/sec) 

Horizontal  Turbulence 

U  U  (CO) 

g  g 

2a2L      1          9        0 
-  -TTU  u     j                    ft"2  rad/sec 

U           Li     0 

°   i+C^)2 

O 

a 
u 

=  5  ft/sec 

L 
u 

=  600  ft 
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TABLE  III 


MODEL  PERFORMANCE  PREDICTIONS 


F 

Baselin 
1.0 

e  Display 
0.5 

0.25 

Qu 
1.0 

ickened  I 
0.5 

)isplay 
0.25 

(ft/sec) 

.881 

1.50 

2,69 

.666 

1.23 

2.45 

a 
w 

(ft/sec) 

2.94 

3,75 

5.11 

2.24 

3.02 

4.39 

a3 
(rad/sec) 

.0109 

.0147 

.0204 

.00829 

.0118 

.0178 

(rad) 

,0104 

.0158 

.0259 

.00795 

.0128 

.0231 

(ft) 

1,81 

2.55 

4.01 

1,34 

2.01 

3.36 

6B 
(ft) 

.0276 

.0300 

,0335 

.0256 

.0277 

.0311 

6C 
(ft) 

.0513 

.0569 

.0675 

.0460 

.0511 

.0611 
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TABLE  IV 


SIMULATION  PERFORMANCE  DATA 


Baseline 

Qui< 

2kened 

Subject 

1 

2      I 

1      1 

2 

a 
u 

5.58t 

2.60 

3.91 

1,86 

(ft/sec) 

1.55tf 

.627 

.653 

.713 

a0 

.0505 

.0269 

.0169 

.0121 

(rad/sec) 

.00482 

.00291 

,00159 

.00134 

ae 

,0471 

.0292 

.0182 

,0197 

(rad) 

,00462 

,00403 

.00614 

.00389 

°h 

9.31 

5.51 

6.88 

4,79 

(ft) 

1,04 

,423 

.949 

.902 

X 

,0534 

.0259 

.0247 

.0127 

(ft) 

.00522 

.00335 

,00269 

.00125 

% 

.0529 

.0373 

.0401 

.0292 

(ft) 

.00642 

.00712 

.00514 

.00501 

t  mean 

. 

tt  standar 

d  deviation 
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r° — 

TELETYPE 


XDS- 9300 

DIGITAL 
COMPUTER 


CI-  5000 

ANALOG 
COMPUTER 


ADAGE  AGT-IO 
DPR-  2 
GRAPHICS 
PROCESSOR 


LINE 
PRINTER 


Figure    1.       Lines   of    communication 
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Figure    2.       Apparatus 
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Units 
20  degrees/inch 
0.5  knots/degree 


Display  Element         Situation 

1.  horizon  bar  pitch  attitude 

2.  airspeed  indicaton    airspeed  error 

3.  aircraft  symbol        stationary 

4.  glideslope  bar        glideslope  deviation   50  feet/inch 

5.  glideslope  indicators  stationary  +  or  -  50  feet 


Figure  3.   Cockpit  display 
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APPENDIX  A  -  DIGITAL  PROGRAM 


INTEGER  IGD<6),FRAME(20),GSLP(12),ACREF<  20) ,HORIZ(25) 

DIMENSION  ITDt60),ITEXT(12),UG( 1000) ,WG( 1000) 

NAMEL 1ST    MGDE,E16,E17,E21,E23,E24,E25 

FRAMEd  )     =    IHEAD(0f6) 
* 

DO    110    1=2,12 

READ    (5,270)     X,Y, IDM 
110    FRAMEd)     =     IPACK(X,Y,  IDM) 

GSLP(l)     =     IHEAD(0,7) 
* 

DO    120    1=2,11 

READ    (5,270)    X,  Y, IDM 
120    GSLP(I)     =     IPACK(X,Y, IDM) 

ACREF(l)     =    IHEAD(0,9) 
* 

DO    130     1=2,19 

READ    (5,270)    X,Y,IDM 
130    ACREF(I)    =    IPACK(XfY, IDM) 

HORIZ(l)     =    IHEAD(0,8) 

DO    140     1=2,5 
READ    ( 5,270)     X, Y, IDM 
140    HORIZ(I)    =    IPACK(XfYtlDM) 
Al=4.472 
A2=3.536 
A3=2.236 
A4=2.738 
A5=2.236 
CMU1=0.07 
OMU2=0.21 
OMU3=0.35 
0MU4=0. 768 
0MU5=1. 765 
0MW1=0.42 
0MW2=1 .26 
0MW3=2. 1 
OMW4=4.608 
OMW5=10.05 
UDR=9.95 
* 

*  TAKE    CARE    OF     PERIPHERAL    EQUIPMENT  * 

OUTPUT( 102) ' ENGAGE    PATCHBOARDS , SELECT    INPUT    CONTROL, EX 
1ECUTE    GATED, AND    SET    SENSE     SWITCH' 

OUTPUT ( 102) 'SELECT    SETTINGS    FOR    POTS    016,017,021,023,0 
124,025' 

OUTPUT( 102) 'AS    E16=    ,E17=    ,E21=     ,E23=    ,E24=    ,E25=    ■ 

INPUT (101) 
* 

CALL  SETPOT  ( 4H? 000 ,. 0091 , 4HP00 1 , .0053 , 4HP002 ,. 0044 ,4H 
1P004,.5680,4HP00  5,.5895,4HP006,.3000,4HP007,.506  7,4HPO 
210,.97S0,4HP011,. 1 6  39 , 4HP01 2 , . 06  00 , 4HP3 1 3 , . 0 1 78 , 4HP0 14 
3,.4732,4HP015,. 15  00, 4HP 0  20, .25  3  4,4HP02  2 , .6667 ,4HP0  27 , . 
47880, 4HP030, .642  0, 4HP031, .3552, 4HP032 , . 6 800 ,4HP033 , . 15 
569,4HP0  37, .0  166, 4HP 04 2,. 3  099 ,4HP045, . 140  0,4HP05  0, .7  725 
6,4HP0  51 , .0242, 4HP052,. 0815, 4HP0 53 ,.1173,4HP054,.7077,4 
7HP056, . 4000, 4HP0 57,. 3500, 4HP025 , E25, 4HP02 1 , E21,4HP023, 
8E23,4HP024,E24,4HP016,Ei6,4HP017,E17) 
IDEV    =    1 

145    IF(SENSESWITCHd)  )     155,150 

150    1F( SENStSWITCHi 2) )     132,145 

152    IDEV=2 

155    CONTINUE 
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CALL    DGINIT    (  IDE  V  ,  IGO  ,  6  ,  I  ER  ) 
CALL    DTINIT    ( IOEV , ITD,60» IER) 

160    OUTPUT( 102) ' SELECT    DISPLAY    MODE,     MODE=i     OR    2« 
INPUT! 101) 

ENCODE     (48,280, ITEXT) 

CALL    TEXTO    (  IDEV  ,  IT  EXT , 12  , 1  0,  1  ,  2 , 2 , I ER ) 
ENCODE     (48,290, (TEXT) 

CALL    TEXTO    (  I DE V , I  TEXT , 1 2 , 1 2 , 1 , 2 , 2, 1 ER) 
ENCODE     (48,300, ITEXT) 

CALL    TEXTO    ( I DEV ,  I  TEXT , 12 , 1 4 , 1 , 2 , 2, 1 ER) 
ENCODE     (48,310, I  TEXT) 

CALL    TEXTO    (  IDE  V ,  I  TEXT , i 2 , 1 6 , 1 , 2 , 2 , 1 ER) 
ENCODE     (48,320, ITEXT) 

CALL    TEXTO    (  IDEV ,  I  TEXT , 12 , 20 , 1 , 2 , 2 , 1 ER ) 
ENCODE     (46,330, I  TEXT) 

CALL    TEXTO    ( I DEV , I  TEXT , 12 , 22 , 1 , 2  ,  2 , 1 ER) 
ENCODE     (48,340, I  TEXT) 

CALL    TEXTO    (  IDEV , I  TEXT , 12 , 24, 1 , 2 , 2 , 1 ER) 
T  =  0.0 

DO    163    1=1,1000 

SU1=0MU1*T 

SU2=OMU2*T 

SU3=0MU3*T 

SU4=0MU4*T 

SU5=0MU5*T 

SW1=0MW1*T 

SW2=OMW2*T 

SW3=0MW3*T 

SW4=0MW4*T 

SW5=OMH5*T 

UG(I)«( l«/25.00)*< A1*SINCSU1)+A2*SIN(SU2)+A3*SIN(SU3)+ 
1A4*SIN( SU4)+A5*SIN(SJ5) ) 

WG( I)=( l./2  5.00)*( A1*SIN( SW1) +A2*SIN( SW2 )+A3*SIN< SW3 )+ 
1A4*SIN( SW4)+A5*SIN(SW5) ) 

T=T+(1.0/UDR) 
163    CONTINUE 
165    IF    (TEST(  D.GT.O)     GO    TO    165 

CALL    DTINIT     1 1  DE  V  ,  I  TD,  60,  IER  ) 

AB    =    0.0 
* 

DO  170  1=6,10 
HQRIZU  )  =  0 
J  =  1+7 

FRAME(J)  =  0 
170  CONTINUE 

IF  (MODE.GT.l)  GO  TO  200 

*  MODE    1,     BASIC    DISPLAY  * 

*********************  5,***  *********************  *************** 

CALL  WRITECLOCK  (0.) 

CALL  RESET  (500) 

CALL  COMPUTE 

CALL  STARTCLOCK 
* 

DO  180  1=1,7200 

CALL    READCLOCK     (V) 

IFCV.GT. 5400.0)    GO    TO    190 

AB    =    AB+1.0 

UG1=UG( I) 

WG1=WG( I) 

CALL    DAC{ L.UGli 2fWGl ) 

CALL    ADK    (OfASi 1 , P , 2 , GS , 3 , Q ,4, GSD j 

P    =    ,2865*P 

GS    =    0.4*GS 

Y  =    .1-P 

H0RIZ(2)     =    IPACK(-.6,Y,0) 
H0RIZ(3)    =    IPACK(~0.125,Y,1) 
H0RIZ(4)    =    IPACKC0.125fYtO) 
H0R1Z(5)     =     IPAC*( .6,Y, 1) 

Y  =    -0.2-GS 
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GSLP(IO)  =  IPACK(-0.57,Y,0) 
GSLP(ll)  =  IPACK (-0.37,Y, 1 » 
GSLP(12)  =  IPACK(-0.57fY,ll 
ANGL  =  0.7029*AS 
SAN  =  SIN(ANGL) 
CAN  =  COS( ANGL) 

Y  =  0.1*SAN+0.1 
Yl  =  Y+0.015 

Y2  =  Y-0.015 

Y3  =  0. 125*SAN*0.1 

X  =  -0. 1*CAN-0.125 

XI  =  X+0.015*SAN 
X2  =  X-C.015-SAN 

X3  "  -0.125*CAN-0.125 
ACREF(2)  =  IPACK(X,Y,0) 
ACREF<3)  =  IPACK (XI, Yl,l) 
ACPEF(4)  =  IPACK(X3,Y3, 1) 
ACRFF<5)  =  IPACK(X2,Y2,1) 
ACREF(6)  =  IPACK(XtYfl) 

X  =  -X 

XI  =  -XI 
X2  =  -X2 
X3  =  -X3 

ACREF(15)  =  IPACK(X,Y,1) 

ACREF(16)  =  IPACK(X2,Y2, 1) 

ACREF(17)  =  IPACK(X3,Y3,i) 

ACREF(18)  =  IPACKCXlt Yltll 

ACREF(19)  =  IPACK(X,Y,1) 

CALL  GRAPHO  ( I DE V , FRAME, 20, 1 , I ER ) 

CALL  GRAPHO  ( I DE V , GSL P , 12 , 2 , I ER ) 

CALL  GRAPHO  ( I DEV , ACkEF , 23, 3, I ER ) 

CALL  GRAPHO  ( IDE  V , HOR 1 Z ,25 ,4, IER ) 
180  CONTINUE 
190  CALL  HOLD  (500) 

GO  TO  260 
* 

*      MODE  2,  QUICKENED  DISPLAY  * 

200  CALL  WRITECLOCK  (0.) 
CALL  RESET  (500) 
CALL  COMPUTE 
CALL  STARTCLOCK 
* 

DO    210    1=1,7200 

CALL    READCLOCK    (  V) 

IF(V.GT. 5400.0)     GO    TO    220 

AB    =    AB+L.O 

UG1=UG( I) 

WG1=WG( H 

CALL    DAC(  1,UG1,2,WG1) 

CALL    ADK    (0,AS, 1 , P, 2 , GS ,3 ,Q ,4 , GSD) 

P    =    . 2665=M  P+2.760*.4*Q) 

GS    =    0.4*(GS+.2*GSD) 

Y  =    .1-P 

HORIZU)    =    IPACK(-.6,Y,0) 
H0RIZ(3)     =     IPACK(-0.125,Y,1 ) 
H0RIZ(4)     =    IPACK(0.125,Y,0) 
H0RIZ(5)     =    1PACK( .6,Y,1) 

Y  =    —0     ?  —  G  S 

GSLP(IO)     =    IPACK(-0.57,Y,0) 
GSLP(ll)     =     IPACK(-0.37,Y,1) 
GSLP(12)     =     IPACK(-0.57,Y,1) 
ANGL    =    0.7029*AS 
SAN    =    SIN( ANGL) 
CAN    =    COS(ANGL) 

Y  =    0.1*SAN*0.1 
Yl    =    Y+0.015 

Y  2    =    Y-  0 ,  3 1  5 

Y3    =    0.  12  5*SAN*0. 1 

X  =    -0. i-CAN-0.l^5 

XI  =    X«-0.015*SAN 
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X2    =    X-0.015*SAN 
X3    =    -0.125*CAN-0.125 
ACREF(2)    =     IPACK(X,Y,0) 
ACREFC3)    =    IPACK(X1,Y1,1) 
ACREF(4)    =    IPACMX3,Y3,  1) 
ACREF(5)    =     IPACK(X2,Y2,  1) 
ACREF(6)     =    IPACK(X,Y,1) 

X  =    -X 

XI  =  -XI 
X2  =  -X2 
X3    =    — X3 

ACREF(15)     =    IPACK(X,Y,1) 

ACREFU6)    =    IPACK(X2,Y2,  I) 

ACREF(17)     =     IPACK(X3,Y3,1 ) 

ACREF(18)    =    IPACK(XlfYltl) 

ACPEF( 19)     =    IPACK(X,Y,  1) 

CALL  GRAPHO  ( I DE V , FRAME , 20 , 1 , I E R ) 

CALL  GRAPHC  ( I DE V , GSL P , 1 2 , 2 , I ER ) 

CALL  GRAPHO  ( I DEV , ACR EF , 20 , 3, I E R ) 

CALL  GRAPHO  ( IDE V ,HOR I Z ,25 , 4, I E R ) 
210  CONTINUE 
220  CALL  HOLD  (500) 
* 

260  CONTINUE 

*  THIS  SECTION!  TAKES  INTEGRATED  SQUARE  VALUES  FROM  THE  * 

*  ANALOG  AND  GIVES  ROOT  MEAN  SQUARE  PERFORMANCE  VALUES.- 

CALL  DGINIT  ( IDEV , IGD t 6, I ER) 

CALL  ADK(l6f SASi 17 , SGS , 20 , SDC , 2 1 , SDB, 22 , THD , 23 , THE  ) 

CALL  RESET  (500) 

SAS  =  SAS/(.104*E25) 

SGS  =  SGS/(.012*E21) 
SDC  =  SDC/(480.*E23) 
SDB  =  SDB/(120.*E24) 
THD  =  THD/(360.0*E17) 
THE  =  THE/(2250.0*E16) 
SAS  =  ABS(SAS) 
SGS  =  ABS(SGS) 
SDC  =  ABS(SDC) 
SDB  =  ABS( SDB) 
THD  =  ABS( THD) 
THE  =  ABS(THE) 
RMAS  =  SQRT(SAS) 
RMGS  =  SQRT(SGS) 
RMDC  =  SQRT(SDC) 
RMDB  =  SQRT(SDB) 
RMTD  =  SQRT(THD) 
RMTH  =  SQRT( THE) 
UDR  =  AB/120.0 
WRITE  (6,350)  MODE 
WRITE  (6,360)  UDR 

WRITE  (6,370)  ^MAS , RMGS , RMDB , RMDC , RMTD, RMTH 
GO  TO  160 
* 

270  FORMAT  (2F10.4,Il) 

280  FORMAT  (•  THIS  IS  A  HELICOPTER  TRACKING 

1    ') 
290  FORMAT  (•       PROBLEM  REQUIRING  THE  ADJUSTMENT  OF 

1    •) 
300  FORMAT  (»       POWER  AND  PITCH  ATTITUDE  TO  MAINTAIN 

1    ■  ) 
310  FORMAT  (»       60  KTS  ON  A  6  DEGREE  GLIDESLOPE. 

1    M 
320  FORMAT  (•  WHEN  READY  TO  BEGIN,  PRESS  THE 

1    ') 
330  FORMAT  (»       RED  BUTTON  ON  THE  COLLECTIVE  LEVER. 

1    •) 
340  FORMAT  (»       THE  TASK  WILL  LAST  FOR  90  SECONDS. 

1    •) 
350  FORMAT  (  ■  1 «  -,  •  DI  S  PLAY  MODE   ■  *  1 1/  ) 
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360  FORMAT  {'O't'UPQATE  RATE  AVERAGED  SFB^,1  TIMES  PER  S 

1EC0NDV) 
370  FORMAT  (,0','RMS  AIRSPEED 


1C 

V'O 

2 

•RMS 

3 

•RMS 

4 

•RMS 

5 

•RMS 

6 

•PMS 

* 


T  'RMS  HEIGHT 

CYCLIC 

COLLECTIVE 

VERTICAL  TURBULENCE 

PITCH  RATE 

PITCH  ANGLE 


»F8.5i 
»F8.5f 
»F8.5, 
tF8.5, 
,F8.5, 


•,F8.5,'     FT/SE 
F8.5,     FT'/'O1, 
FT'/'0», 
FT'/« 0' , 
FT/SEC'/'O' , 
RAD/SEC'/'O1 , 
RAD  •/) 


END 
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PROGRAM  VARIABLES 

AB  Incremental  number  of  display  updates  during  simulation 

ACREF  Vector  of  coordinates  defining  aircraft  reference. 

ANGL  Angle  defining  airspeed  deviation. 

AS  Airspeed  deviation. 

CAN  Cosine  of  given  angle. 

E16  Setting  for  potentiometer  016. 

E17  Setting  for  potentiometer  017. 

E21  Setting  for  potentiometer  021. 

E23  Setting  for  potentiometer  023. 

E24  Setting  for  potentiometer  024. 

E25  Setting  for  potentiometer  025. 

FRAME  Vector  of  coordinates  defining  display  frame. 

GS  Glideslope  (height)  deviation. 

GSD  Glideslope  deviation  rate. 

GSLP  Vector  of  coordinates  defining  glideslope  reference. 

HORIZ  Vector  of  coordinates  defining  horizon  reference. 

IDEV  Graphics  terminal  in  use,  1  or  2 . 

IDM  Graphics  instruction,  draw  or  move. 

IER  Error  type  code , 

IGD  Graphics  instruction,  size  of  graphics  input. 

ITD  Graphics  instruction,  size  of  text  input. 

ITEXT  Vector  of  text  input. 

MODE  Display  mode,  1  or  2. 

P  Pitch  angle. 
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Q  Pitch  rate. 

RMAS  Root  mean  square  airspeed  deviation. 

RMDB  Root  mean  square  cyclic  movement. 

RMDC  Root  mean  square  collective  movement. 

RMGS  Root  mean  square  glideslope  deviation. 

RMTD  Root  mean  square  pitch  rate. 

RMTE  Root  mean  square  pitch  angle, 

SAN  Sine  of  given  angle, 

SAS  Mean  square  airspeed  deviation. 

SDB  Mean  square  cyclic  movement. 

SDC  Mean  square  collective  movement, 

SGS  Mean  square  glideslope  deviation. 

T  Time  increment, 

THD  Mean  square  pitch  rate. 

THE  Mean  square  pitch  angle. 

UDR  Update  rate  for  the  90  second  period. 

UG  Vector  of  horizontal  turbulence  magnitudes, 

UG1  Specific  horizontal  turbulence  value. 

WG  Vector  of  vertical  turbulence  magnitudes. 

WG1  Specific  vertical  turbulence  value. 

V  Internal  clock  count,  60  cycles  per  second. 
X  Running  variable,  abscissa. 

Y  Running  variable,  ordinate. 

A1,A2,...A5  Sine  wave  magnitudes. 

0MU1,0MU2, , . ,0MU5    Horizontal  turbulence  frequencies 

0MW1,0MW2, , . .0MW5    Vertical  turbulence  frequencies. 
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SU1 ,SU2, , . , SU5         Sine  values  for  horizontal  turbulence 
SW1,SW2 , . . , SW5         Sine  values  for  vertical  turbulence. 
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APPENDIX  B  -  ANALOG  PROGRAM 


W 


-Hi 


[i] 
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